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We measure the local contact forces at both the top and bottom boundaries of three-dimensional 
face-centered-cubic and hexagonal-close-packed granular crystals in response to an external force 
applied to a small area at the top surface. Depending on the crystal structure, we find markedly 
different results which can be understood in terms of force balance considerations in the specific 
geometry of the crystal. Small amounts of disorder are found to create additional structure at both 
the top and bottom surfaces. 
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I. INTRODUCTION 

In a static granular bead pack all of the forces on a 
particle, due to the contacts from its neighbors and due 
to gravity, must be in perfect balance; there can be no 
net force on any particle. Moreover, in the limit of hard 
particles, no particle distortion is allowed leading several 
investigators to question whether elasto-plasticity theory 
can be used in this hard-sphere limit to describe the me- 
chanical properties of an amorphous granular material 
[|l|-Q . Others suggest that elastic theory is capable of pro- 
ducing the observed experimental results [|]|]] . The pres- 
ence of disorder makes an exact calculation of the contact 
forces in such a material impossible. On the other hand, 
it should be possible to calculate the forces exactly in a 
perfect crystal |7j]. However, because macroscopic grains 
do not spontaneously crystallize, it has been difficult to 
create large defect-free granular crystals ]^-pj| on which 
to perform an experiment to see how the forces are trans- 
mitted. In addition, there will always be small imperfec- 
tions or asperities in the particles themselves and slight 
variations in the diameters of the grains. Although we 
cannot completely avoid disorder of this latter kind, we 
have found an efficient method of constructing essentially 
perfect crystals with few, if any, defects in the particle 
positions With such crystals, one would hope that 
the presence of small disorder from the variation of the 
individual particle shapes and sizes could be treated as a 
perturbation about the perfect crystalline response. 

Previously, there have been experimental studies of the 
statistical properties of the local contact forces in granu- 
lar packs. In particular, the distribution of normal forces, 
P(F), between the particles and the container walls was 
measured in a variety of situations [|l2|-[l|]. In most of 
those studies, a piston pressed down on the entire top 
surface of the granular material. The surprising result 
was that the entire form of P(F) (a large value near 
F = 0, a small peak [or plateau] with a maximum near 
the average force and an exponential tail at large forces) 
was robust and did not depend on whether the granu- 
lar medium was amorphous or crystalline or whether the 



particles were smooth or rough 0] . The question is thus 
raised as to what aspects of the force propagation in a 
granular material are affected by preparation and the un- 
derlying crystal or amorphous structure. 

Several experiments have studied the response (i.e., 
the Green's function) to a localized external force at the 
top surface of a granular medium. The different stud- 
ies found somewhat contradictory results. DaSilva and 
Rajchenbach [[l6| studied a system consisting of a spe- 
cial two-dimensional packing of photoelastic rectangular 
bricks. They measured the contact forces between par- 
ticles as a function of position and found that the con- 
tact forces had a maximum centered below the position 
of the applied perturbation and that on ether side the 
forces decreased in magnitude. The width of this peak 
in the forces grew as the square root of the depth below 
the surface. This result is consistent with models pre- 
dicting a diffusive propagation of the forces such as the 
q- model proposed by Coppersmith et al. |l7],[l^]. Other 
studies have focused on the response to a localized ex- 
ternal force in amorphous packings of spheres. Rey del- 
let and Clement found that the response at the bottom 
of a three-dimensional pack had a central maximum of 
approximately Lorenztian shape with a width that grew 
lineraly with the depth of the pack j^J . Similar behavior 
was found in two-dimensional samples studied by Geng 
et al. |2(J . This linear increase of the width with depth is 
consistent with elastic theories ]|||. The work by Geng 
et al. also suggests a strong effect due the spatial or- 
dering of the particles resulting in a central minimum of 
the response forces in ordered packings. Thus there is 
a need to study the propagation of forces from a local- 
ized perturbation in a nearly perfectly ordered crystalline 
array of particles in order to investigate how the underly- 
ing order (or disorder) of the lattice influences the force 
propagation. 

Bouchaud et al. have proposed a model in which 
forces propagate in straight lines until encountering a 
defect where they are split, and they have shown that 
this model leads to a central force minimum for small 
depths and to elastic-like behavior for large pack depths 



1 



The straight-line propagation of forces (resulting 
from the hyperbolic model of Cates et al. ^| ) and the 
splitting of forces can be viewed in a geometrical frame- 
work based upon the contact orientations of neighboring 
beads. From this viewpoint crystalline packings allow for 
the control of these contact orientations and the resulting 
splitting. In particular, the hexagonal-close-packed (hep) 
crystal structure is such that forces are required to split 
at nearly every layer as they travel downward through 
the pack. By looking at the response of an hep crystal to 
a localized force it is possible to explore the consequences 
of multiple splittings in comparison with the predictions 
of the forcechain splitting model. 

In this paper we present an experimental investiga- 
tion of the response to a localized external force in three 
dimensional nearly perfect hep and face-centered-cubic 
(fee) crystals. We measure the resulting contact forces 
at both the top and bottom surfaces of the crystal, and 
investigate how the pattern of forces depends on the crys- 
tal structure and the depth of the pack. We determine 
that in the regime of highly ordered granular packs, the 
crystal structure plays a dominant role in determining 
the average pattern of response forces at the bottom sur- 
face. These patterns can be viewed as a consequence of 
force balance in the specific geometry of the crystal. We 
also consider the effect of small amounts of disorder in 
these crystals and show that they modify the pattern of 
forces on the bottom surface in a nontrivial manner and 
scatter forces back upwards through the pack to the top 
surface. 



II. EXPERIMENTAL METHODS 

We created large nearly perfect hep and fee granu- 
lar crystals of 3.06 ± 0.04 mm diameter soda lime glass 
spheres. These crystals of approximately 53 bead diam- 
eters on a side were contained within an acrylic cylin- 
der with a triangular cross-section designed to be com- 
mensurate with the hep crystal structure. The crystals 
were created as individual layers of triangular order in the 
manner described in reference (T^|. In this way the po- 
sition of each individual bead was viewed and corrected 
if neccessary as the crystal was constructed, so that the 
resulting crystal would have very few defects (estimated 
to be less than ten beads significantly out of place for a 
crystal of approximately 50000 beads). Furthermore, the 
stacking order of the planes was controlled in order to 
achieve three dimensional fee and hep crystals. 

After building such a crystal we applied an external 
force to a localized region at the center of the top of the 
crystal either by using a slowly applied steady force of 
approximately 100 pounds or by using a quick impulse 
force. When using the slowly applied force it was nec- 
essary to apply the force over a region approximately 
six beads in diameter. Attempts to apply the force to 
a smaller region resulted in significant bead movement. 
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FIG. 1. (Color) A comparison of the pattern of calibrated 
forces in response to a slowly applied force to the pattern 
of intensities in response to a quick impulse force, (a) The 
pattern of calibrated forces at the bottom surface of a 19 
layer fee crystal in response to a slowly applied force over an 
area six beads in diameter at the top surface, (b) The pattern 
of intensities of carbon marks in response to a quick impluse 
force applied to the same crystal structure over the same area. 
Both pictures represent an average over ten experimental runs 
and an average over crystal symmetries. 

When using the quick impulse force it was possible to 
apply the force to a smaller area of approximately two 
beads in diameter. For the purpose of comparison ex- 
periments were also done with the quick impulse force 
applied to the large six bead diameter region. 

We used a carbon paper technique to measure the 
forces in response to the external forces ]T^ , p!3| , p!7| , p3[ . A 
piece of carbon paper and a piece of white paper (24|] 
were placed on the top and bottom of the crystal. As the 
force was applied individual beads pressed into the car- 
bon paper and left marks on the white paper. The size 
and darkness of the resulting marks was related to the 
normal force on the corresponding beads. When using 
the slowly applied force, we calibrated the force versus 
darkness response as explained in [p^| . The positions of 
the carbon marks were matched to a triangular lattice 
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and the forces for each lattice point were averaged over 
approximately ten experimental runs. It is important to 
note that individual runs vary significantly, but that the 
average over many runs remains highly reproducible. In 
order to improve statistics, the locally measured forces 
were also averaged over symmetries of the crystal (one 
reflection and two rotations). 

Note that the analysis technique must be carried out 
carefully so as not to introduce error from the averaging 
over many files, and crystal symmetries. If the centers 
of the images are not determined accurately, then lattice 
positions from file to file will not match up exactly and 
will result in a small degree of broadening of any features. 
We estimate this broadening to be significantly less than 
one bead diameter for each of the regions of large force. 

When using the quick impulse, we calibrate the mea- 
sured forces by comparing them with those obtained 
with an identical set-up using the slowly applied force 
described above. We compared the patterns of carbon 
marks left by the two methods and found no qualitative 
difference. Figure |l](a) shows the pattern of calibrated 
forces at the bottom of a 19 layer fee crystal in response 
to a slowly applied force at the top surface over a region 
six beads in diameter. The calibration maps points of 
zero intensity to a small but nonzero force which results 
in a nonzero background. Figure 0(b) shows the pattern 
of intensities of carbon marks at the bottom of a similar 
crystal in response to a quick impulse force applied to the 
same area. The color scale of intensity has been shifted 
by a factor of ten and a minimum value used to match 
the previous pattern of forces. Using this shift in color 
scale intensity we map the two different measurements 
onto one another and obtain a calibration for the quick 
impulse. 

III. RESULTS 

In amorphous packings it is difficult to apply a local- 
ized force large enough to leave carbon marks without 
also having significant movement of the grains. As such 
we were only able to explore the response of amorphous 
packings to the quick impulse force applied to a large 
six bead diameter region. In this case we observed the 
pattern of intensities at the bottom of three-dimensional 
amorphous packs of various depths and found a broad 
central maximum. Our data on the width of this central 
maximum are consistent with a square-root growth with 
depth as seen in a two-dimensional packing of bricks by 
DaSilva and Rajchenbach jlfj, but are also within error 
of a linear growth as seen by Reydellet and Clement |HJ 
in three-dimensional packings. 
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FIG. 2. (Color) Intensity patterns at the bottom surface of 
fee crystals (a) 9 and (b) 19 layers in depth in response to a 
quick impulse force applied to a two bead diameter region at 
the top surface. 

A. FCC crystals 

When using highly ordered packings bead movement 
is much less of a problem, and it is possible to use a 
slowly applied force to a large area or a quick impulse 
to a smaller area. For fee crystals we found a charac- 
teristic response pattern at the bottom surface which is 
qualitatively independent of type of force applied. Fig- 
ure [| shows the patterns of intensities at the bottom of 
fee crystals, 9 layers and 19 layers in depth, in response 
to a quick impulse force applied to an area two beads in 
diameter at the center of the top surface of the crystal 
averaged over approximately ten experimental runs and 
averaged over crystal symmetries. We find three regions 
of large force, arranged in a triangular pattern. The size 
of this triangle increases with the depth of the crystal. 
The size of the three regions of large force are compa- 
rable to the size of the region over which the force was 
applied. In addition to the three regions of large force, 
there appear fainter but still significant regions of force 
along the lines connecting these three regions. 

Figure ||(a,b) shows cross sections of the intensity ver- 
sus position across the triangular pattern at the bottom 
surface of 9 and 19 layer fee crystals. As depicted in 
the inset of figure |3|(a) the cross-section proceeds across 
the edge of the triangular pattern, corresponding to a 
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FIG. 3. Intensity as a function of position at the bottom of 
FCC crystals (a) 9 and (b) 19 layers in depth in response to a 
quick impulse force applied to an area two beads in diameter 
at the top of the pack. The cross- sect ion was taken through 
the triangular pattern of intensity as shown in the inset, (c) 
The position of main intensity peaks as a function of the depth 
of the pack in response to (circles) a quick impulse force over 
a small area, (triangles) a quick impulse force over a large 
area, and (squares) a slowly applied force over a large area. 
Error bars represent the full width at half the maximum value 
of the peaks. The black line indicates a slope of 35.3 degrees 
with respect to the vertical which corresponds to an angle 
specific to the geometry of the fee crystal. The inset shows 
the position of the smaller peaks corresponding to the edge 
of the triangular pattern with a line representing an angle of 
19.5 degrees. 
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FIG. 4. (Color) Intensity patterns at the bottom of hep 
crystals (a) 9 and (b) 25 layers in depth in response to a 
quick impulse force over an area two beads in diameter at the 
top surface of the crystal. 

small peak, through the center of the pattern (position 
= 0) with small intensities, and across the region of large 
force, corresponding to a large peak. The position of the 
large peaks as a function of the depth of the crystal are 
indicated in figure |](c), with the error bars corresponding 
to the full width at half maximum of the peaks, which 
can in some cases be seen to be significantly skewed. Ad- 
ditional data points are included in part c corresponding 
to experimental runs with a quick impulse applied to a 
large area and to runs with a slowly applied force over 
a large area. The sloped line indicates an angle of 35.3 
degrees which corresponds to the angle down the edge of 
a tetrahedron. The regions of large force are found to lie 
just inside of this angle as explained in section The 
inset of figure |](c) shows the position of the small peaks 
corresponding to the edge of the triangular patterns as 
a function of the depth of the crystal packing. The solid 
line shows an angle of 19.5 degrees with respect to the 
vertical, which is corresponds to the angle down the face 
of a tetrahedron. 



B. HCP crystals 

In contrast to the triangular patterns found for fee 
crystals, hep crystals show rings of large force at the 
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FIG. 5. Radial distribution of intensity at the bottom of 
hep crystals (a) 9 and (b) 25 layers in depth in response to a 
quick impulse force applied to a small area at the top surface. 
The position of the peaks as a function of the depth of the 
pack are shown in (c) as circles with error bars representing 
the full width at half of the maximum value of the peak. 
As the interior of the ring is nearly uniform no lower bound 
is shown for the 25 layer crystal. The triangular data points 
indicate similar peaks in response to a quick impulse force over 
a large area, and the square data points correspond to a slowly 
applied force over a large area. The solid line corresponds to 
an angle of 27.4 degrees with respect to the vertical which is 
the average of two angles specific to an hep crystal structure 
as explained in the text. 

bottom surface in response to a localized force at the 
top surface. Figure || shows the intensity patterns at the 
bottom surfaces of hep crystals 9 and 25 layers in depth 
in response to a quick impulse force applied to an area 
two beads in diameter at the top surface averaged over 
approximately ten experimental runs and averaged over 
crystal symmetries. As the depth of the pack increases 
the radius of the ring increases, and interior fills in. Fig- 
ure ||(a,b) shows the radial distribution of intensity for 
hep crystals of 9 and 25 layers in depth in response to a 
quick impulse force over a small area. Figure g(c) shows 
the position of the peaks in intensity as a function of 
the depth of the pack with the error bars corresponding 
to the width of the peak at half of the maximum value. 
For hep crystals 25 layers in depth the intensity is nearly 
uniform and thus no lower bound is shown. Additional 



data points are included corresponding to experimental 
runs with a quick impulse applied to a large area and to 
runs with a slowly applied force over a large area. The 
sloped line indicates an angle of 27.4 degrees which is 
the average of the angles down the edge and the face of a 
tetrahedron. The peaks in intensity and force are found 
to lie just inside of this angle as explained in section [Tv| . 

C. Top surface 

In addition to measuring the response at the bottom 
of the crystal packs, we also used the carbon paper tech- 
nique to measure the response at the top of the pack. 
Figure ^ shows surface plots for the average intensity 
of the carbon dots as a function of position at the top 
surface of the crystal in response to a quick impulse ap- 
plied to an area two beads in diameter. Both fee and 
hep crystals show a characteristic six-fold spoke pattern, 
which does not depend significantly on the depth of the 
crystal. The magnitude of the average intensity in these 
spokes is approximately ten times less than the intensity 
in the regions of large force on the bottom surface, but 
they are still significantly stronger than the surrounding 
background. The spokes extend all of the way to the 
edge of the analysis region (15 bead diameters in radius) 
without a significant decrease in intensity. 

IV. GEOMETRICAL INTERPRETATION 

A. FCC geometry 

Many aspects of the patterns of force or intensity seen 
in figures [l|, ||, [|, and || can be explained through force 
balance in the specific geometry of the crystals. Our 
fee crystals were oriented with the 111 plane horizontal, 
so that they consisted of a stack of triangularly ordered 
planes. This is equivalent to a typical cannon ball piling. 
In an fee crystal the stacking order from plane to plane 
is such that every third layer lies on top of the first. This 
results in straight lines of contacts between beads running 
diagonally down through the pack. 

Figure |?] shows the geometry of a piece of a perfect fee 
crystal. If a downward force is applied to the bead at the 
top of this pyramidal pile, then the force is transmitted to 
the three beads below it in the second layer. Essentially, 
the bead at the top of the pyramid is supported by a 
tripod, with three lines of beads stretching downward 
through the pack. Only the beads in this tripod will 
feel a force in response to the force on the top bead. 
All other beads inside and outside of the pyramid are 
shielded from that force. Thus, at the bottom of the 
crystal there are only three beads which feel any force. 
These beads are arranged in a triangular pattern with 
the size of the triangle proportional to the depth of the 
pack. The legs of the tripod in the fee crystal are at an 
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FIG. 6. (Color) Intensity patterns at the top surface of 
(a) a 3 layer hep crystal, (b) a 9 layer hep crystal, (c) a 9 
layer fee crystal, and (d) a 19 layer fee crystal in response to 
a quick impulse force applied to a small area in the center 
of these images. Each image is averaged over approximately 
ten experimental runs and averaged over crystal symmetries. 
The six-fold spoke pattern is present in both types of crystal 
packings and does not depend on the depth of the packing. 



angle of 35.3 degrees with respect to the vertical. Given 
this angle the position of the beads which feel the applied 
force is determined. The solid line in figure ||(c) shows 
the predicted angle of 35.3 degrees. 

B. disorder 

If there is a small amount of disorder in the crystal, 
then the positions of the beads within the tripod legs 
will not be perfectly aligned. Thus, in order to balance 
the forces vectorally, there must be some force scattered 
off of the tripod. For a small degree of disorder, the 
majority of the force will still be transmitted along the 
tripod legs, with a small amount scattered to one or more 
additional neighboring beads. Figure [t|(c) represents an 
fee crystal with one point of disorder in the second layer. 
Most of the force is transmitted along the tripod legs, 
but a small amount is scattered along the blue line. If 
there is no further disorder, then this scattered force will 
propagate along its own line of beads all of the way to 
the bottom of the pack, and arrive at the bottom surface 
along one of the lines connecting the three ends of the tri- 
pod legs, as shown in figure @(d). Using this as a model 
for the effect of one point of disorder we find that, when 
considering a small degree of disorder throughout the sys- 
tem, the majority of the force is transmitted along the 
tripod legs, with a small amount scattered down the sides 
of the pyramid and arriving at the bottom surface along 
the edges of the triangle formed by the tripod legs. Fig- 
ure ^(e) shows the first order effect of scattering within 
an fee crystal. The red lines represent strong force, and 
the blue lines represent weaker forces resulting from first 
order scattering. The pattern of forces at the bottom of 
the crystal, shown in ^|(f) is seen in the images in figures 
and | 

Disorder can also explain the appearance of the regions 
of large force at angles slightly smaller than 35.3 degrees 
as seen in figure ||(c) as scattered forces arriving preferen- 
tially inside the triangular pattern will skew the regions 
of large force towards the inside of the pyramid. To first 
order, no force exists outside of the pyramid. 

If at points of disorder the force scatters to beads in 
plane rather than to the plane below, then without addi- 
tional scattering the scattered force propagates straight 
outward to the sides of the container. For nonclose- 
packed systems, in which beads in-plane are not in con- 
tact, it is also possible for the forces to scatter upwards 
through the pack. This would produce a spoke-like struc- 
ture, but the length of the spokes would be expected to 
depend on the depth of the pack which is not seen exper- 
imentally. 

The second-order effect of disorder on the transmission 
of forces is depicted in figure ||. Figure ||(a) shows a piece 
of a three layer fee crystal. The red lines represent the 
strong forces along the tripod discussed earlier. The blue 
lines are the result of first order scattering with both in- 



6 




FIG. 7. (Color) Geometric interpretation of the response to 
a point force in an fee crystal, (a) In a perfect fee crystal the 
force at the top of the pyramid is balanced by a tripod, (b) 
At the bottom layer only three beads support force, (c) One 
point of disorder along a tripod leg results in a scattering of 
a small amount of force off of the tripod leg. In this case the 
force is shown being scattered downward (d) This scattered 
force propagates downward to the bottom surface to one of the 
lines connecting the regions of large force, (e) Averaging over 
all possible points of single downward scattering events, the 
scattered forces are contained on the faces of the pyramidal 
structure, (f) The scattered forces arrive at the bottom layer 
arranged along the lines connecting the non-scattered forces. 




FIG. 8. (Color) (a) Schematic diagram of an fee crystal 
showing second order scattering. The red lines represent large 
unscattered forces. Blue lines correspond to first order scat- 
tering, both downward and in-plane, and green lines represent 
some of the second order scattering. Also shown are the pat- 
terns at the (b) top and (c) bottom surfaces of a 5 layer fee 
crystal due to second order scattering. 

plane and downward scattering shown. The green lines 
represent some of the scattering from the scattered forces, 
or second order scattering. Of interest is the scattering 
from the in-plane forces. These forces can scatter up- 
wards, downwards, or in-plane through the pack. Figure 
||(b) shows the forces at the top surface of a five layer 
deep fee crystal, where second order scattering produces 
a six fold spoke pattern extending radially outward from 
the center after a one-bead gap. These spokes extend 
outward to the edge of the system and do not depend on 
the depth of the crystal. Similar experimental spoke pat- 
terns are shown in figure ||. Figure ||(c) shows the pattern 
of forces at the bottom surface of the crystal where the 
effect of second order scattering is to fill in the triangle 
with very weak forces and extend the edges of the trian- 
gle beyond the vertices. For small amounts of disorder 
the magnitude of the forces from second-order scattering 
is much smaller than those previously discussed. Nev- 
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FIG. 9. Schematic diagrams of a perfect hep crystal, (a) 
Forces cannot be transfered from the top to the third layer 
without having negative forces, shown as dashed lines. When 
the weight of the particles is neglected this pack is unstable 
without additional support from the sides, (b) Illustration of 
force being transfered to in-plane neighbors 



ertheless figure [|(a) shows faint forces along the lines of 
the triangle edges beyond the vertices. 



C. HCP Geometry 

The hep crystals that we studied are also oriented so 
as to consist of stacked planes of triangular order. Unlike 
the fee crystals the stacking of the planes is such that ev- 
ery other plane lies directly above the first. The top two 
layers of an hep crystal and fee crystal are identical, and 
the force transmission will be the same. However the dif- 
ference in position of the beads in the subsequent layers 
causes the force transmission to be markedly different. 
In an hep crystal there are not straight lines of beads 
downward through the pack, so even in a perfect crystal 
the force cannot be transferred to a single bead and still 
preserve vector force balance. Furthermore, the nature 



of the geometry of the hep crystal is such that the force 
on a bead in the second layer received from the top bead 
cannot be transferred to the three beads below. This 
would require a negative force (see figure ^)) or adhesion, 
which is not present in dry granular systems. The result 
is that a pack such as that shown in figure §(a) is not 
stable to an external force when there is no friction and 
the weight of the individual grains is neglected. If only 
the weight of the beads is considered then a perfect hep 
packing is stable down to five layers after which it devel- 
ops the same instability. In order to keep the beads from 
moving there must be additional support. Inside a larger 
crystal this instability demands that forces are transmit- 
ted to beads in-plane for close-packed crystals and that 
forces are transmitted upwards for nonclose-packed crys- 
tals. If forces are allowed to be transmitted to neighbors 
in-plane then figure ^|(b) shows the propagation of forces 
for a small piece of a perfect hep crystal in response to 
a downwards force at the top. In-plane forces will con- 
tinue directly to the walls without additional scattering 
in a perfect crystal. First order scattering will result in a 
spoke like pattern at the top surface similar to that seen 
in the fee crystals. The expected pattern at the bottom 
surface is a ring of maximum force, whose radius grows 
with additional layers of depth at angles of 19.5 degrees 
and 35.3 degrees for alternating layers. The solid black 
line in figure || shows the expected radius of the ring as a 
function of depth based on an average angle of 27.4 de- 
grees. As the depth of the pack increases the interior of 
the ring also fills in. This process occurs even in a perfect 
hep crystal. Disorder has little effect on the pattern of 
forces on the bottom surface, but does result in a spoke 
pattern on the top surface, as is seen experimentally in 
figure pi 



V. CONCLUSIONS 

We have observed that in the highly ordered regime 
the response to a localized force is greatly dependent 
on the structure of the beads within a granular pack. 
We found distinctly different responses from fee crystals 
and hep crystals, both of which behave differently than 
amorphous packs. Amorphous packs are found to re- 
spond with a large central peak, fee crystals respond 
with three regions of large force at the bottom surface, 
while hep crystals show rings of maximum force. These 
structures can be explained by force balance in the ge- 
ometry of the crystal. These findings are consistent with 
the two-dimensional results of Geng et al. |20[ , |25|l . 

The texture of a granular pack is crucially important 
for determining the average pattern of stress transmis- 
sion, as has been suggested by Rajchenbach |26|. The fee 
crystal structure, without disorder, leads to a tripod of 
force chains that is the dominant feature of the force pro- 
files. Small amounts of disorder create secondary struc- 
tures at the bottom surface and are responsible for a six- 
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fold spoke pattern of forces at the top surface. The hep 
crystal structure does not allow for straight line propaga- 
tion of forces downward through the pack. Even without 
disorder the forces must split at every layer. 

The splitting of forces within an hep pack can be 
viewed in analogy to the forcechain splitting model of 
Bouchaud et al In a perfect hep crystal the force 

splitting is driven by the orientations of neighboring 
grains and is thus nearly everywhere uniform, in con- 
strast to the forcechain splitting model in which ran- 
domly placed defects drive forces to split into random 
angles. Nevertheless, the average response of the sys- 
tem to multiple splittings is found to have similarities 
across these two systems. The forcechain splitting model 
has been shown to lead to a central minimum of force 
at the bottom surface in response to a localized force 
at the top for small packing heights, which is consistent 
with our observations on hep crystals of small heights. In 
two-dimensional simulations of forcechain splitting with 
a small density of defects the response pattern made a 
transition towards elastic-like behavior with a broad cen- 
tral peak for depths larger than approximately six mean 
free paths of the force We find that for three- 

dimensional hep crystals larger than approximately ten 
layers in depth the central minimum of force fills in lead- 
ing to a broad central peak response. This experimen- 
tal result in three dimensions with splitting generated 
uniformly throughout the pack by contact orientations 
agrees qualitatively with the predictions of the forcechain 
splitting model simulations in two dimensions where ran- 
domly placed defects generate the splitting. 

Finally, it is important to note that previous exper- 
iments have shown that the crystal structure does not 
affect P(F), the probability distribution of interparticle 
contact forces fll2| ]. We conclude that the fluctuations of 
forces within a single pack are dominated by microscopic 
details of the granular pack and are not influenced by 
the arrangement of the beads. Nevertheless, the average 
response to a localized force is greatly influenced by this 
structure. 
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